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Abstract 

An inclusive search for anomalous production of two prompt, isolated leptons with the same electric 
charge is presented. The search is performed in a data sample corresponding to 4.7 fb -1 of integrated 
luminosity collected in 2011 at ^ = 7 TeV with the ATLAS detector at the LHC. Pairs of leptons 
(e ± e ± , e ± fi ± , and A i± ^ ± ) with large transverse momentum are selected, and the dilepton invariant 
mass distribution is examined for any deviation from the Standard Model expectation. No excess is 
found, and upper limits on the production cross section of like-sign lepton pairs from physics processes 
beyond the Standard Model are placed as a function of the dilepton invariant mass within a fiducial 
region close to the experimental selection criteria. The 95% confidence level upper limits on the cross 
section of anomalous e ± e ± , e ± /x ± , or /i ± / u ± production range between 1 .7 fb and 64 fb depending on 
the dilepton mass and flavour combination. 
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1 Introduction 

Events containing two prompt leptons with large transverse momentum (px) an d equal 
electric charge are rarely produced in the Standard Model (SM) but occur with an enhanced 
rate in many models of new physics. For instance, left-right symmetric models [1-4], Higgs 
triplet models [5-7], the little Higgs model [8], fourth-family quarks [9], supersymmetry [10], 
universal extra dimensions [11], and the neutrino mass model of refs. [12-14] may produce 
final states with two like-sign leptons. In the analysis described here, pairs of isolated, 
high-px leptons are selected, and the invariant mass of the dilepton system (e^e^, e =t /u ± , 
is examined for the inclusive final state and separately for positively- and negatively- 
charged pairs. 

The ATLAS Collaboration has previously reported inclusive searches for new physics in 
the like-sign dilepton final state in a data sample corresponding to an integrated luminosity 
of 34 pb" 1 [15] and in like-sign muon pairs with 1.6 fb _1 [16]. No significant deviation from 
SM expectations was observed, and fiducial production cross-section limits as well as limits 
on several specific models of physics beyond SM were derived. The CDF Collaboration has 
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performed similar inclusive searches [17, 18] without observing any evidence for new physics. 
Furthermore, the ATLAS and CMS Collaborations have performed several searches for like- 
sign leptons produced in association with jets or missing transverse momentum where no 
evidence of non-SM physics was observed [19-26]. 

This article is organised as follows. A brief description of the ATLAS detector is given 
in section 2. The data and simulation samples, the event selection, and the background 
determination are explained in sections 3, 4, and 5, respectively. The systematic uncer- 
tainties on the background estimate (including theoretical uncertainties on the production 
cross sections) and signal acceptances are summarised in section 6. In section 7 the number 
of observed lepton pairs in data is compared to the background estimate, and in section 8 
these results are used to derive upper limits on the fiducial cross section for like-sign dilep- 
ton production in a kinematic region closely related to the experimental event selection. 

2 The ATLAS detector 

The ATLAS detector [27] consists of an inner tracking system, calorimeters, and a muon 
spectrometer. The inner detector (ID), directly surrounding the interaction point, is com- 
posed of a silicon pixel detector, a silicon microstrip detector, and a transition radiation 
tracker, all immersed in a 2 T axial magnetic field. It covers the pseudorapidity 1 range 
1 77 1 < 2.5 and is enclosed by a calorimeter system consisting of electromagnetic and hadronic 
sections. The electromagnetic part is a lead/liquid-argon sampling calorimeter, divided into 
a barrel (\r]\ < 1.475) and two end-cap sections (1.375 < \r/\ < 3.2). The barrel (|r/| < 0.8) 
and extended barrel (0.8 < \rj\ < 1.7) hadronic calorimeter sections consist of iron and scin- 
tillator tiles, while the end-cap (1.5 < |r/| < 3.2) and forward (3.1 < \r/\ < 4.9) calorimeters 
are composed of copper or tungsten, and liquid-argon. 

The calorimeter system is surrounded by a large muon spectrometer (MS) built with 
air-core toroids. This spectrometer is equipped with precision tracking chambers (com- 
posed of monitored drift tubes and cathode strip chambers) to provide precise position 
measurements in the bending plane in the range |r/| < 2.7. In addition, resistive plate 
chambers and thin gap chambers with a fast response time are used primarily to trigger 
muons in the rapidity ranges |r/| < 1.05 and 1.05 < |r/| < 2.4, respectively. The resistive 
plate chambers and thin gap chambers also provide position measurements in the non- 
bending plane, which are used for the pattern recognition and the track reconstruction. 

The ATLAS trigger system has a hardware-based Level- 1 trigger followed by a software- 
based high-level trigger [28]. The Level-1 muon trigger searches for hit coincidences between 
different muon trigger detector layers inside geometrical windows that define the muon 
transverse momentum and provide a rough estimate of its position. It selects high-px muons 
in the pseudorapidity range |r/| < 2.4. The Level-1 electron trigger selects local energy 

1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in 
the centre of the detector and the z-axis along the beam line. The i-axis points from the IP to the centre 
of the LHC ring, and the y-axis points upward. Cylindrical coordinates (r, <fi) are used in the transverse 
(x, y) plane, <fi being the azimuthal angle around the beam line. The pseudorapidity is defined in terms of 
the polar angle 6 as r\ — — lntan(6>/2). 
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clusters of cells in the electromagnetic section of the calorimeter. For low-energy electron 
clusters, low activity in the hadronic calorimeter nearby in the rj-cf) plane is required. The 
high-level trigger selection is based on similar reconstruction algorithms as those used 
offline. 

3 Data sample and Monte Carlo simulation 

This analysis is carried out using a data sample corresponding to an integrated luminosity 
of 4.7 ± 0.2 fb -1 of pp collisions, recorded in 2011 at a centre-of-mass energy of 7 TeV. In 
this dataset, the average number of interactions per beam crossing ranges from about six 
in the first half of the year to about fifteen at the end of 2011. 

The data were selected using single-muon and single-electron triggers with px thresh- 
olds of 10 GeV and 16 GeV at Level-1, respectively. In the high-level trigger, a muon with 
px > 18 GeV is required, while for electrons, the px threshold is 20 GeV in the early 2011 
data and 22 GeV in the later part of the year. To ensure no efficiency loss for electrons 
with very high-px, a trigger with a px threshold of 45 GeV is also used which has no 
requirement on the hadronic calorimeter energy deposits near the electron in the rj-cfr plane 
at Level-1. 

Monte Carlo (MC) simulation is used to estimate some of the background contributions 
and to determine the selection efficiency and acceptance for possible new physics signals. 
The dominant SM processes that contribute to prompt like-sign dilepton production are 
WZ and ZZ, with smaller contributions from like-sign W pair production (W^V^ 1 * 1 ) and 
production of a W or Z boson in association with a top quark pair (tiW, and tiZ). These 
are all estimated using MC simulation. For processes with a Z boson, the contribution 
from 7* — > £ + £~ due to internal or external bremsstrahlung of final-state quarks or leptons 
is also simulated for m(£ + £~) > 0.1 GeV. Sherpa [29] is used to generate WZ and ZZ 
events, while W ± W ± , ttW, and ttZ production is modelled using MadGraph [30] for the 
matrix element and Pythia [31] for the parton shower and fragmentation. The W^W^ 
sample includes the W^W^W^ process. 

The normalisation of the WZ and ZZ MC samples is based on cross sections de- 
termined at next-to-leading order (NLO) with MCFM [32]. The cross sections times 
branching ratios for W^Z — > l^vt^i^ and ZZ — > i^F^l^l^ after requiring two charged 
leptons (electrons, muons, or taus) with the same electric charge and with px > 20 GeV 
and \rj\ < 2.5, are 372 fb and 91 fb, respectively. 

For tiW and ttZ production, the higher-order corrections are calculated in ref. [33] 
and refs. [34-36], respectively. The full higher-order corrections for VF^VF 1 * 1 production 
have not been calculated. However, for parts of the process, the NLO QCD corrections 
have been shown to be small [37, 38]. Based on this, no correction is applied to the LO 
cross section. 

Opposite-sign dilepton events due to Z/'y*, ti, and W^W^ production constitute a 
background if the charge of one of the leptons is misidentified. The Z/'y* process is gen- 
erated with Pythia, and the cross section is calculated at next-to-next-to-leading order 
(NNLO) using PHOZPR [39]. The ratio of this cross section to the leading-order cross 
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section is used to determine a mass dependent QCD K-factor which is applied to the result 
of the leading-order simulation [40]. Higher-order electroweak corrections (beyond the pho- 
ton radiation included in the simulation) are calculated using Horace [41, 42], yielding 
an electroweak if-factor due to virtual heavy gauge-boson loops. The production of tt is 
modelled using MC@NLO [43], with Herwig [44] used for parton showering and hadro- 
nisation. The normalisation is obtained from approximate NNLO QCD calculations using 
Hathor [45]. The production of W ± VF =F is generated using Herwig, and the cross section 
is normalised to the NLO value calculated with MCFM. The production of Z7 — > iij and 
Wj — > ivy constitutes a background if the photon converts to an e + e~ pair. Alpgen [46] 
is used to generate Wj production, and the cross section is normalised to be consistent 
with the ATLAS Wj cross-section measurement for pt(t) > 15 GeV [47]. The Zj process 
is implicitly accounted for in simulation of the Z/j* process. 

In addition, a variety of new physics signals are simulated in order to study the effi- 
ciency and acceptance of the selection criteria. Pair production of doubly-charged Higgs 
bosons (pp — > i? ±± if =F=F ) via a virtual Z/j* exchange is generated using Pythia for H ±i: 
mass values between 50 GeV and 1000 GeV [48]. Production of a right-handed W bo- 
son (Wr) decaying to a charged lepton and a right-handed neutrino (Nr) [49, 50] and 
pair production of heavy down- type fourth generation quarks (1^4) decaying to tW are 
generated using Pythia. Like-sign top-quark pair production can occur in models with 
flavour-changing neutral currents, e.g. via a t-channel exchange of a Z' boson with utZ' 
coupling. Since the left-handed coupling is highly constrained by B® — B® mixing [51], only 
right-handed top quarks (tn) are considered. A sample is generated with the Protos [52] 
generator, based on an effective four-fermion operator uu — > tt corresponding to Z' masses 
S> 1 TeV [53]. The parton shower and hadronisation are performed with Pythia. 

Parton distribution functions (PDFs) taken from CTEQ6L1 [54] are used for the Mad- 
Graph and Alpgen samples, and the MRST2007 LO** [55] PDF set is used for the 
Pythia and Herwig samples. For the tt MC@NLO sample, CTEQ6.6 [56] PDFs are 
used, while for the generation of diboson samples with Sherpa the CTEQ10 [57] parame- 
terisation is used. Uncertainties on the quoted cross sections are discussed in section 6. 

The detector response to the generated events is simulated with the ATLAS simulation 
framework [58] using Geant4 [59], and the events are reconstructed with the same software 
used to process the data. Simulated minimum bias interactions generated with Pythia 
are overlaid on the hard scatter events to closely emulate the multiple pp interactions 
present in the current and in adjacent bunch crossings (pileup) present in the data. The 
simulated response is corrected for the small differences in efficiencies, momentum scales, 
and momentum resolutions observed between data and simulation. 

4 Event selection 

Events are selected if they contain at least two leptons (ee, e/i, ///i) of the same electric 
charge with > 20 GeV. If the higher-px lepton is an electron, it is required to have 
Pt > 25 GeV. At least one muon with pr > 20 GeV or one electron with px > 25 GeV is 
required to match an object that passed the relevant high-level trigger. These px thresholds 
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are chosen such that they are on the plateau of the trigger efficiency. In addition to the 
lepton selection, events must have a primary vertex reconstructed with at least three tracks 
with pt > 0.4 GeV. If more than one interaction vertex is reconstructed in an event, the 
one with the highest X^tj summed over the tracks associated with the vertex, is chosen 
as the primary vertex. For Z bosons decaying to electrons or muons, the efficiency of the 
primary vertex requirements is close to 100%. 

Electrons are identified as compact showers in the electromagnetic sections of the 
calorimeters that are matched to a reconstructed track in the ID using the tight criteria 
described in ref. [60]. These criteria include requirements on the transverse shower shapes, 
the geometrical match between the track and shower, the number and type of hits in 
the ID, and they reject electrons associated with reconstructed photon conversion vertices. 
Electrons must have \rj\ < 2.47, excluding the transition region between the barrel and end- 
cap calorimeters (1.37 < < 1.52). The electron tracks are refitted with the Gaussian 
Sum Filter algorithm to account for radiative energy losses (bremsstrahlung) in the detector 
material [61]. 

Muon candidates are formed from tracks reconstructed in the ID combined with tracks 
reconstructed in the MS [62]. To reduce the charge mismeasurement rate the independent 
charge measurements from these two detectors are required to agree. Muons must have 
\<n\ < 2.5. 

The ID tracks of electron and muon candidates are required to be consistent with 
originating from the primary interaction vertex. The requirements are |e?o| < 1-0 mm, 
[zosin^l < 1.0 mm, and \do/a(do)\ < 3, where do (zq) is the transverse (longitudinal) 
distance between the primary vertex and the point of closest approach of the track and 
cr(do) is the estimated uncertainty on the do measurement, typically 15 /an. To ensure 
precise impact-parameter measurements, a hit in the innermost pixel detector layer is 
required if the track traverses an active detector component. This requirement also reduces 
background due to electrons from photon conversions. 

Both electrons and muons are required to be isolated from other activity in the event. 
The isolation energy, sum of transverse energies 2 in calorimeter cells (in- 

cluding electromagnetic and hadronic sections) in a cone AR = \J (At/) 2 + (Ac/>) 2 < AR\ so 
around the lepton direction. This quantity is corrected for the energy of the electron object 
as well as energy deposits from pileup interactions. The track isolation, p™ ncARis ° ( j s analo- 
gously defined as the scalar sum of the transverse momentum of tracks with p^ > 1 GeV in 
a cone of AR < AR{ so around the lepton direction, excluding the lepton track. For muons, 
pCone0.4^ T ^ ^ q Qg anc j ^concO.4 < 4 GeV + 0.02 x pr(n) are required. For electrons, the 
requirements are pf nc0 - 3 /p T (e) < 0.1 and £™ne0.2 < 3 GeV + (pr(e) - 20 GeV) x 0.037. 
These requirements are chosen to maintain good efficiency for prompt leptons over a large 
Pt range while rejecting a large fraction of non-prompt leptons, particularly at low pt- 

In addition to the isolation criteria, selected leptons are required to be well sepa- 
rated from jets to suppress leptons from hadronic decays. For this purpose, jet candidates 

2 Transverse energy is the projection of an energy deposition in the calorimeter in the transverse plane, 
where the direction is given by the line joining the energy deposition and the nominal interaction point. 
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are reconstructed from topological clusters in the calorimeter [63] using the anti-fcj algo- 
rithm [64] with a radius parameter of 0.4. Their energies are corrected for calorimeter 
non-compensation, energy loss in upstream material, and other instrumental effects. Fur- 
thermore, quality criteria are applied to remove reconstructed jets not arising from hard- 
scattering interactions [65]. Jets are required to have pr > 25 GeV, \rj\ < 2.8, and 75% of 
the momentum of the tracks within the jet must originate from the primary vertex. Any 
jet within a distance AR = 0.2 of a candidate electron is removed from the list of jets to 
avoid counting the electron as a jet. Electrons and muons are then required to be separated 
by AR > 0.4 from any jet with p T > 25 GeV + 0.05 x pr(£). 

Lepton pairs are selected if they contain two leptons with the same electric charge 
passing the above selection requirements with invariant mass m(£ ± £ ± ) > 15 GeV. Fur- 
thermore, for e ± e ± pairs, the mass range 70-110 GeV is vetoed due to large backgrounds 
from opposite-sign electron pairs produced by Z decays, where the charge of one lepton 
is misidentified. Any combination of two leptons is considered, allowing more than one 
lepton pair per event to be included. 

5 Background determination 

Backgrounds to this search arise from three principal sources: SM production of prompt, 
like-sign lepton pairs (prompt background); production of opposite-sign lepton pairs where 
the charge of one lepton is misidentified (charge-flip background) or where a photon pro- 
duced in association with the opposite-sign leptons converts into an e + e~ pair; and hadrons 
or leptons from hadronic decays (non-prompt background). 

5.1 Prompt lepton backgrounds 

SM processes resulting in two prompt leptons of the same electric charge are WZ, ZZ, 
W^W^, tiW, and tiZ production. The background due to these processes is determined 
from MC simulation using the samples and cross sections described in section 3. Other SM 
processes are not expected to contribute significantly to the background and are neglected. 

5.2 Charge misidentification and photon conversions 

Charge misidentification can occur for high-momentum tracks when the tracking detector 
is unable to determine the curvature of the track. This misidentification can occur for 
both electrons and muons. In the momentum range of a few hundred GeV relevant for this 
analysis, the background arising from this source is determined to be negligible based on 
MC simulation. The rate of misidentification is studied in data with muons from Z — > \x\i 
decays where the muons are selected based only on the MS (ID) measurements to study 
the ID (MS) mismeasurement. The combined probability to mismeasure the charge in both 
the ID and the MS is found to be consistent with zero and an upper limit ranging from 
5 x 10~ 10 at low to 7 x 10~ 2 for j>t ~ 400 GeV is placed. This upper limit is used to 
determine the systematic uncertainty for muons. 

Another source of charge misidentification occurs for electrons when a high-momentum 
photon is radiated and converts into an e + e~ pair. Electron candidates are rejected if 
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they are associated with conversion vertices [60]. However, for asymmetric conversions, 
sometimes only one of the tracks is reconstructed and its charge may be different from 
the charge of the original electron that radiated the photon. The probability for electron 
charge misidentification to occur, either by this mechanism or by track mismeasurement, 
is measured as the fraction of like-sign e ± e ± pairs with 80 GeV< m(ee) < 100 GeV as a 
function of the electron rj following the procedure explained in ref. [66]. The simulation 
is found to overestimate this probability by about 15%. An 77-dependent scaling factor is 
applied to the simulated rate to account for this overestimate, and the simulation is then 
used to predict the backgrounds from Z/j*, ti, and M /± M /=F production. In the central 
(forward) region the probability increases from about 0.01% (0.6%) at px ~ 50 GeV to 
0.3% (4%) at pt ~ 300 GeV. An uncertainty on this misidentification probability of ±(10- 
20)% (depending on rj) is derived by comparing different methods to determine this factor 
as described in ref. [66] and by considering the statistical uncertainty on each method. 

Production of W"f — > ivy and Zj — > £ + £~^f can lead to like-sign lepton pairs if the pho- 
ton converts. This background is determined using the MC samples described in section 3. 
Since this background is closely related to the charge misidentification for electrons de- 
scribed above, the MC-based estimate is also scaled down by 15% and the same systematic 
uncertainty of ±(10-20)% is applied. 

5.3 Non-prompt lepton backgrounds 

The non-prompt background includes lepton pairs where one or both of the leptons result 
from hadronic decay or misidentification. Processes contributing to this background include 
W+jet, Z+jet, multi-jet (including bb and cc), and ti production. 

The non-prompt background is determined directly from data. For both electrons and 
muons, the determination relies on measuring a factor / that is the ratio of the number 
of selected leptons satisfying the analysis selection criteria [N$) to the number of leptons 
failing these selection criteria but passing a less stringent set of requirements, referred to 
as anti-selected leptons {NjCj- 

f= Ns_ 
1 N A 

The factor / is determined as a function of pt and 77 in a data sample dominated by non- 
prompt leptons. Any contamination by prompt leptons in the numerator and denominator 
is subtracted using MC simulation. The details of the anti-selected definitions and the 
regions used for measuring / depend on the lepton flavour and are described in more detail 
below. 

The primary sources of non-prompt muons are semi-leptonic decays of b- and c-hadrons. 
For the anti-selection, the isolation criterion described in section 4 is inverted, but a looser 
requirement of p™ ne0 a /pt(/^) < 1-0 is placed. The measurement of / is done in a sample 
dominated by non-prompt muons with \do\/a(do) > 5 and |g?o| < 10 mm. These muons 
are selected from dimuon events with pt(m) > 10 GeV, m(fifi) > 15 GeV, and where 
Z-boson candidates have been vetoed. Simulated muons from b- and c-hadron decays 
are used to determine the expected difference between / for muons with \do\/a(do) > 5 
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(used in the data-driven determination of /) and muons with \do\/a(do) < 3 (used in 
the event selection). In simulation, muons with \do\/a(do) < 3 are more isolated, so a 
correction factor of 1.34 ± 0.34 is applied to /. The resulting value of / is about 0.10 at 
Pt = 20 GeV and increases to about 0.25 at = 100 GeV. One source of uncertainty in 
this procedure is that the fraction of non-prompt muons from pion and kaon decays may 
be larger for muons passing the analysis-level impact parameter criteria than for those 
used in the measurement of /. The uncertainty due to the contribution of light-flavour 
decays is assessed by exploiting the difference in the track momentum measurement in 
the ID and the MS, which is expected to be large for kaons or pions decaying in the 
ID or calorimeters. The difference between / for light-flavour decays and heavy- flavour 
decays is combined with the fraction of non-prompt muons from light-flavour decays to 
give an uncertainty of ±15% on /. The total systematic uncertainty on / is the quadratic 
sum of the statistical error, the uncertainty on the prompt background subtraction, the 
correction for the isolation-dependence on the impact parameter significance criteria, and 
the uncertainty in the contribution of pion and kaon decays. It is about ±37% at low pt 
and ±100% for p T > 100 GeV. 

For electrons, the main non-prompt backgrounds arise from heavy-flavour decays, 
charged pions that shower early in the calorimeter, and neutral pions decaying to two 
photons where one of the photons converts to an e + e~ pair. The anti-selected electrons 
must have either \do\/a(do) > 3 or fail the medium electron identification criteria while 
passing the loose electron criteria, as defined in ref. [60]. The inversion of the impact pa- 
rameter criteria enhances the heavy-flavour background, while the inversion of the medium 
identification criteria enhances the pion backgrounds. A sample dominated by non-prompt 
electrons is selected by triggering on a high-pT electron candidate and requiring that the 
event contains a jet with pt > 20 GeV. Contamination of this sample by prompt electrons 
is reduced by requiring that there be only one electron candidate and requiring the trans- 
verse mass 3 formed by the electron and the missing transverse momentum (E^f lss ) to be 
below 40 GeV. The value of / is found to be about 0.18 at Ei = 20 GeV and decreases to 
about 0.10 at E? = 100 GeV. A systematic uncertainty is determined accounting for the 
statistical uncertainty, the prompt background subtraction, the pt requirement of the jet 
in the event, and any trigger bias. The uncertainty is about ±10% at low E<i and increases 
to ±100% for £t > 300 GeV. Alternate anti-selection definitions are used to assess the 
overlap between the non-prompt and charge misidentification backgrounds, as well as the 
contributions of light- and heavy-flavour jets to the non-prompt background. For the for- 
mer, a correction is made to the non-prompt background, amounting to ±(11-25)%, and 
the full correction is taken as a systematic uncertainty. For the latter, an uncertainty of 
±(4— 84)% is found, depending on the invariant mass threshold considered. 

A background prediction is derived from / using dilepton pairs where one or both 
leptons are anti-selected but pass all other event selection criteria. The non-prompt back- 

3 The transverse mass is defined as y / 2E™ BB pT(e)(l — cos A<f>(e, i5™ lss )), where _E™ 1SS is the missing trans- 
verse momentum and A0(e, _E™ 1SS ) is the difference between the azimuthai angles of the electron and E™ lss . 
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ground prediction is then given by 

Na+s N,s+a Na+a 

^non-prompt = /(PTl,T/l) + / (PT2 , ?72 ) ~ ^ /(PTl, »7l)/ (PT2,?72), (5.1) 

i=l i=l i=l 

where the sums are over the number of pairs, iV^+s (iVs+,4), where the subleading (leading) 
lepton satisfies the selection criteria and the other lepton satisfies the anti-selection, or over 
the number of pairs where both leptons are anti-selected (Na+a)- The prompt contribu- 
tions to Na+s, Ns+a, and Na+a are subtracted based on MC predictions. The variables 
Pti anci rji (pt2 and 772) refer to the kinematic properties of the leading (subleading) lepton. 

The background estimates are verified in several control regions that are designed to 
probe specific sources of background. The understanding of prompt leptons is tested in a 
sample requiring two leptons of opposite charge. This sample is dominated by the Z/j* 
process in all three final states, and the data agree with the background prediction to better 
than 5%. The kinematics of the electron charge misidentification background are verified 
in pairs of like-sign electrons with 80 GeV < m(e ± e ± ) < 100 GeV. The non-prompt 
background is tested by inverting the isolation criteria but requiring a looser isolation 
cut value. The heavy-flavour background is tested by inverting the \do\/a(do) criterion. 
Furthermore, for electrons, the light-flavour background is tested by inverting some of the 
shower shape criteria. The data agree with the background prediction in the various control 
regions within the systematic uncertainties. 

6 Systematic uncertainties 

Several systematic effects can change the signal acceptance and the background estimate. 
Uncertainties on the event selection efficiencies and on the luminosity affect the predicted 
yield of signal events as well as those backgrounds that are estimated from MC simula- 
tion. Uncertainties on the trigger efficiencies, lepton identification efficiencies, and lepton 
momentum scales are determined using W, Z, and J/ip decays [67-69]. The uncertainties 
on the lepton identification efficiencies result in an uncertainty on the number of like-sign 
pairs of ±(3-4)%. Effects resulting from uncertainties in the muon momentum scale and 
resolution are negligible. For electrons the energy scale is known to ±1%, and its impact 
on the signal acceptance is below 0.1%, while the background estimate is affected by up 
to ±3.5%. The uncertainties due to the trigger efficiencies are < 1% in all channels. The 
uncertainty on the integrated luminosity is ±3.9% [70, 71]. 

The uncertainties on the production cross sections of the SM processes affect the pre- 
dicted yield of the prompt lepton background. For WZ and ZZ production, an uncertainty 
of ±10% is estimated due to higher-order QCD corrections by varying the renormalisation 
and factorisation scales by a factor of two. The uncertainty due to the parton distribution 
functions (PDFs) is evaluated using the eigenvectors provided by the CTEQ10 set [57] of 
PDFs, following the prescription given in ref. [72]. The difference between the central cross- 
section values obtained using this PDF set and that obtained with the MRST2008NLO 
PDFs [55] is also added in quadrature. This procedure gives a conservative estimate of 
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the PDF uncertainty on the cross section. The resulting uncertainty on the cross section 
is ±7%. For tiW, tiZ and W^W^ production a normalisation uncertainty of ±50% is 
assigned [34, 37]. 

The production of Z/"/*, tt, and W also constitutes a background in the ee 

and e/i channels due to charge misidentification. The PDF uncertainties for Z/7* and ti 
are computed using the MSTW2008 NNLO PDF sets [73] and the renormalisation and 
factorisation scales are varied by factors of two to derive the uncertainties due to higher- 
order QCD corrections. For Z/7*, an additional systematic uncertainty is attributed to 
electroweak corrections [40]. The total uncertainties on the Z/j* and it cross sections are 
±7% and +10/-11%, respectively. The W^W* cross-section uncertainty is determined 
using the same scale and PDF variations as for WZ and ZZ, yielding a total uncertainty 
of ±12%. 

The systematic uncertainty on the background from non-prompt leptons is discussed 
in section 5.3. It includes a component due to the systematic uncertainties associated with 
the measurement of the ratio / and a component due to the statistical uncertainty on the 
number of anti-selected leptons. For m(^ ± ^ ± ) > 15 GeV it is ±33% for the ±31% for 
the e/i, and ±28% for the ee final states, and it increases to nearly ±100% at higher masses 
for all final states. 

The probability of assigning the wrong charge to the lepton and its uncertainty is 
discussed in section 5.2. For the /i/x final state, charge misidentification results in an uncer- 
tainty of to 9 events for m(^ ± ^ ± ) > 15 GeV and ±^ 7 events for ?n(£ ± ^ ± ) > 400 GeV. For 
electrons, the uncertainty on the charge- flip background due to uncertainty in the charge 
mismeasurement probability ranges from ±15% to ±23% for the same mass thresholds. 
The uncertainty on the VF7 and Z7 backgrounds from this source is ±(15-18)% depending 
on m{l±l±). 

Statistical uncertainties due to the limited size of the background MC samples are also 
considered. Systematic uncertainties on different physics processes due to a given source 
are assumed to be 100% correlated. 

7 Comparison of data to the background expectation 

The predicted numbers of background pairs are compared to the observed numbers of like- 
sign lepton pairs in table 1. For the e^e^ (e^^) final state, about 30% (50%) of the 
background is from prompt leptons in all mass bins. The rest arises from leptons from 
hadronic decays, charge flips, or photon conversions. For the final state, the prompt 

background constitutes about 83% of the total background for m(// ± /i ± ) > 15 GeV, rising 
to nearly 100% of the total background for m(^ ± /x :l= ) > 400 GeV. The overall uncertainty 
on the background is about ±15% at low mass and ±30% at high mass. 

Table 2 shows the data compared to the background expectation separately for £ + £ + 
and £~ t~ events for the ee, fifj,, and efx final states. The background is higher for the £ + £ + 
final state due to the larger cross section in pp collisions for W + than W~ bosons produced 
in association with 7, Z, or hadrons. 
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Table 1. Expected and observed numbers of pairs of isolated like-sign lcptons for various cuts 
on the dilcpton invariant mass, to(^ ± £ ± ). The uncertainties shown are the quadratic sum of the 
statistical and systematic uncertainties. The prompt background contribution includes the WZ, 
ZZ, W ± W ± , tiW, and ttZ processes. When zero events are predicted, the uncertainty corresponds 
to the 68% confidence level upper limit on the prediction. 



Sample 


Numt 
> 15 GeV 


>er of electron 
> 100 GeV 


pairs with m 
> 200 GeV 


(e ± e ± ) 
> 300 GeV 


> 400 GeV 


Prompt 

Non-prompt 

Lnarge nips and 
conversions 


101 ± 13 
75 ±21 

170 ± 33 


56.3 ±7.2 
28.8 ±8.6 

91 ± 16 


14.8 ±2.0 
5.8 ±2.5 

22.1 ±4.4 


4.3 ±0.7 

5 +0 - 8 

8.0 ±1.7 


1.4 ±0.3 

o o +0 - 2 

3.4 ±0.8 


Sum of backgrounds 


346 ± 44 


176 ± 21 


42.8 ±5.7 


12.8 ±2.1 


4.8 ±0.9 


Data 


329 


171 


38 


10 


3 




Num 
> 15 GeV 


1 r 

ber ot muon t 
> 100 GeV 


>airs with m(f 
> 200 GeV 


i ii ) 
> 300 GeV 


> 400 GeV 


Prompt 
Non-prompt 
Lnarge nips 


205 ± 26 
42 ± 14 

n n+4.9 

O-0-o.o 


90 ±11 
12.1 ±4.6 

n n+2.5 

- -o.o 


21.8 ±2.8 
1.0 ±0.6 

0- -o.o 


5.8 ±0.9 

o o +0 - 3 

n n+1.7 

O-0-o.o 


2.2 ±0.4 

o o +0 - 3 

n n+1.7 

- -o.o 


bum of backgrounds 


247_2g 


102 ± 12 


on c+3.4 
ZZ.O_2 g 


K Q+l-9 

b - 8 -0.9 


2 - 2 -o.4 


Data 




1 1 n 


9Q 


u 


9 




Num 

> 15 GeV 


aer of lepton 
> 100 GeV 


aairs with m( 
> 200 GeV 


> 300 GeV 


> 400 GeV 


Prompt 

Non-prompt 

Charge flips and 
conversions 


346 ± 43 
151 ± 47 

142 ± 28 


157 ± 20 
45 ± 13 

33 ± 7 


36.6 ± 4.7 
9.2 ± 4.1 

10.5 ± 2.8 


10.8 ± 1.5 
2.6 ± 1.1 

2.9 ± 1.2 


3.9 ± 0.6 
1.0 ± 0.6 

2.2 ± 1.1 


Sum of backgrounds 


639 ± 71 


235 ± 25 


56.4 ± 7.0 


16.3 ± 2.3 


7.0 ± 1.4 


Data 


658 


259 


61 


17 


7 



The level of agreement between the data and the background expectation is evaluated 
using 1— CLft [74], defined as the one-sided probability of the background-only hypothesis 
to fluctuate up to at least the number of observed events. Statistical and systematic 
uncertainties and their correlations are fully considered for this calculation. The largest 
upward deviation, observed for > 100 GeV, occurs about 8% of the time in 

background-only pseudo-experiments for this mass bin. 

Figure 1 shows the dilepton invariant mass spectra for the e ± e ± , // ± // ± , and e^ji^ 
final states, and figure 2 shows the pt of the leading leptons. Good agreement with the 
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Table 2. Expected and observed numbers of positively- and negatively-charged lcpton pairs for 
different lower limits on the dilepton invariant mass, m(£ ± ^ ± ). The uncertainties shown are the 
quadratic sum of the statistical and systematic uncertainties. 



Sample 


Number of lepton pairs with m( 
> 15 GeV > 100 GeV > 200 GeV 


t ± i ± ) 
> 300 GeV 


> 400 GeV 






e+e + 


psirs 






Sum of backgrounds 


208 ± 28 


112 ± 14 


28.6 ± 4.0 


8.5 ±1.4 


3.3 ±0.7 


Data 


183 


93 


26 


6 


1 








psirs 






Sum of backgrounds 


138 ± 21 


63.3 ± 8.5 


14.2 ± 2.3 


4.4 ±0.8 


i n a 

1 ^4+ u - 4 
i - O4 -0.3 


Data 


1 AP, 
140 


1 o 


1 9 


4 


2 






1J+IJ-* 


pall o 






Sum of backgrounds 


147 ± 17 


63 7 + H 

U,J 1 —7.6 


14 5+?"i 




1 6 +0 - 9 
i -°-0.3 


Data 


144 


fin 


ID 


4 


2 








pairs 






Sum of backgrounds 


100 ± 12 


38-4+5;° 


8.31H 


1 7+0.9 
L -'-0.3 


6 +0 - 9 
u - D -0.l 


Data 


120 


50 


13 


Z 


U 






eV 


pairs 






Sum of backgrounds 


381 ± 42 


142 ± 15 


33.8 ± 5.3 


9.8 ± 1.5 


4.2 ± 0.9 


Data 


375 


149 


39 


9 


4 








pairs 






Sum of backgrounds 


259 ± 31 


93 ± 10 


22.6 ± 3.0 


6.5 ± 1.3 


2.9 ± 1.0 


Data 


283 


110 


22 


8 


3 



background expectation is observed for all three channels. 

8 Upper limits on the cross section for prompt like-sign dilepton pro- 
duction 

Since the data are consistent with the background expectation, upper limits are placed 
on contributions due to processes from physics beyond the Standard Model. Limits are 
derived independently in each final state and mass bin shown in tables 1 and 2. The 95% 
confidence level (C.L.) upper limit on the number of events, -/V95, is determined using the 
CL S method [74] with a test statistic based on the ratio between the likelihood of the signal 
plus background hypothesis and the likelihood of the background-only hypothesis. The 
likelihoods follow a Poisson distribution for the total number of events in each signal region, 
and are calculated based on the number of expected and observed events. Systematic 
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ATLAS + Data 2011 

jut = 4.7 fb 1 □Non-prompt 
Vs = 7TeV □ Charge flips 
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mfeV) [GeV] 



ATLAS -t- Data 2011 

j"l_dt = 4.7 fb-' □ Non-prompt 
Vs = 7TeV □ P^pt 




300 350 400 
m(nV) [GeV] 



(a) 



(b) 



> 

CD 

CD 



^ J' 



ATLAS 4- Data 2011 

Ldt = 4.7fb- 1 □ Non-prompt 
□ Charge flips 




50 100 150 200 250 300 350 400 
m(eV) [GeV] 



(c) 



Figure 1. Invariant mass distributions for (a) e ± e ± , (b) A t± / i± i arL d (c) e ± /x ± pairs passing the 
full event selection. The data are shown as closed circles. The stacked histograms represent the 
backgrounds composed of pairs of prompt leptons from SM processes, pairs with at least one non- 
prompt lcpton, and for the electron channels, backgrounds arising from charge misidentification and 
photon conversions. Pairs in the ee channel with invariant masses between 70 GeV and 110 GeV 
are excluded because of the large background from charge misidentification in Z — > e ± e =F decays. 
The last bin is an overflow bin. 



uncertainties are incorporated into the likelihoods as nuisance parameters with Gaussian 
probability density functions. For the inclusive f^t^ final states, the upper limit ranges 
from 168 pairs for m(efi) > 15 GeV to 4.8 pairs for m(///x) > 400 GeV. 

The limit on the number of lepton pairs can be translated to an upper limit on the 
cross section measured in a given region of phase space (referred to here as the fiducial 



region), crfg, via 



fid ^95 

£fid x J L-dt 

where is the efficiency for detecting events within the fiducial region and j Cdt is the 
integrated luminosity of 4.7 fb -1 . 

The fiducial region definition is based on MC generator information such that it is 
independent of the ATLAS detector environment. The value of egd generally depends on 
the new physics process, e.g. the number of leptons in the final state passing the kinematic 
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Figure 2. Leading lepton distributions for (a) e ± e ± , (b) ( u ± ^ ± , and (c) e ± ^ ± pairs passing 
the full event selection. The data are shown as closed circles. The stacked histograms represent 
the backgrounds composed of pairs of prompt leptons from SM processes, pairs with at least one 
non-prompt lepton, and for the electron channels, backgrounds arising from charge misidentification 
and photon conversions. The last bin is an overflow bin. 



selection criteria or the number of jets that may affect the lepton isolation. In order to 
minimise this dependence, the definition of the fiducial region is closely related to the 
analysis selection. In the MC generator, electrons or muons are selected as stable particles 
that originate from a W or Z boson, from a r lepton, or from an exotic new particle (e.g. a 
or a right-handed W). The generated electrons and muons are required to satisfy the 
Pt> an d isolation requirements listed in table 3, which mirror the selection requirements 
described in section 4. The generator-level track isolation, p^ neARia °, is defined as the 
scalar sum of all stable charged particles with > 1 GeV in a cone of size Ai?j so around 
the electron or muon, excluding the lepton itself. The requirement on isolation energy 
(£'™ ncARlso ) placed in the electron selection is not emulated in the fiducial region definition 
because the measurement of isolation energy in the calorimeter has a poor resolution. Pairs 
of like-sign leptons must have m(£ ± ^ ± ) > 15 GeV and for e^e^, the mass range 70-110 GeV 
is additionally excluded. 

The fiducial efficiency, egd, is defined as the fraction of lepton pairs passing this se- 
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Table 3. Summary of requirements on generated leptons in the fiducial region. The definition of 
the isolation variable, p^ neARia °, is given in the text. 





Electron requirement 


Muon requirement 


Leading lepton px 


p T > 25 GeV 


p T > 20 GeV 


Sub-leading lepton pt 


p T > 20 GeV 


p T > 20 GeV 


Lepton r] 


\r]\ < 1.37 or 1.52 < |t?| < 2.47 


M < 2.5 


Isolation 


p cone0.3/ pT < q j 


pf ne0A /p T < 0.06 and 

p cone0.4 < 4 Q e y + q Q2 x ^ 



lection that also satisfy the experimental selection criteria described in section 4. It is 
determined for a variety of new physics models: i? ±± bosons with mass between 50 GeV 
and 1000 GeV [48], fourth generation quarks decaying to Wt with masses of 300-500 GeV, 
like-sign top-quark production via a contact interaction [53], and right-handed W bosons 
decaying to £ ± £ ± via a right-handed neutrino with m(Wij) = 800-2500 GeV [49, 50]. These 
models are chosen as they cover a broad range of jet multiplicities and lepton pt spectra. 

For m(£ ± £ :iz ) > 15 GeV in the e ± e ± channel, the fiducial efficiencies range from 43% 
for models with low-px leptons to 65% for models with high-px leptons. The primary 
reason for this dependence is that the electron identification efficiency varies by about 15% 
over the relevant px range [60]. The model dependence introduced by not emulating the 
calorimeter isolation in the definition of the fiducial region is < 1%. For the e^fx^ channel, 
£fid ranges from 55% to 70%, and for the p ± /x ± final state it varies between 59% and 72%. 
For the higher dilepton mass thresholds the efficiencies are slightly larger than for the lower 
mass thresholds. The efficiencies are also derived for £ + £ + and £~£~ pairs separately and 
found to be independent of the charge. For the same new physics models, the fraction 
of events satisfying the experimental selection originating from outside the fiducial region 
ranges from < 1% to about 9%, depending on the final state and the model considered. 

To derive the upper limit on the fiducial cross section, the lowest efficiency values are 
taken for all mass thresholds, i.e. 43% for the e^e 1 * 1 , 55% for the e ± ^ ± , and 59% for the 
/ u ± /Lt ± analysis. The 95% C.L. upper limits on the cross section, Ugg , are calculated using 
equation 8.1 and given in table 4. The limits range from 64 fb to 1.7 fb for the inclusive 
analysis depending on the mass cut and the final state, and are generally within la of the 
expected limits. Upper limits on the £ + £ + and l~ i~ cross sections are also derived, using 
the same efficiencies as for the inclusive limits. 

9 Conclusions 

A search for anomalous production of like-sign lepton pairs has been presented using 
4.7 fb -1 of pp collision data at ^fs = 7 TeV recorded by the ATLAS experiment at the 
LHC. The data are found to agree with the background expectation in e^e^, e^fi^, and 
fi^f-i^ final states both in overall rate and in the kinematic distributions. The data are used 
to constrain new physics contributions to like-sign lepton pairs within a fiducial region of 
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Table 4. Upper limits at 95% C.L. on the fiducial cross section for £ ± £ ± pairs from non-SM physics. 
The expected limits and their la uncertainties are given, as well as the observed limits in data, for 
the ee, e/x, and fifi final state inclusively and separated by charge. 



Mass range 


expected 
e ± 


9 

observed 
e± 


5% C.L. upper limit [ft 
expected | observed 


expected 


observed 


m, > 15 GeV 


461- 


42 


561^ 


64 


24.01^ 


29.8 


m > 100 GeV 


24.1±gJ 


23.4 


23.01^ 


31.2 


12 2 +4 - 5 
iz ' z -3.0 


15.0 


m > 200 GeV 




7.5 


8.41^ 


9.8 


4.31};? 


6.7 


m > 300 GeV 


4 -°-1.3 


3.9 


4 1+1-8 
4 -l-0.9 


4.6 


9 4 +o.9 

z - 4 -0.7 


2.6 


m > 400 GeV 


2 q +L1 

z - y -0.8 


2.4 


°- u -0.8 


3.1 


-i 7+O.6 
l-'-0.5 


1.7 










77i > 15 GeV 


29 1 +10 - 2 
zy - x -8.6 


22.8 


34.9l^ 6 2 


34.1 


15.01^ 


15.2 


m > 100 GeV 


16.lj£§ 


12.0 


15.411;? 


18.0 


8-4±l;l 


7.9 


m > 200 GeV 


7.0111 


6.1 


6-6111 


8.8 


o r+1.6 
°-°-0.7 


4.3 


m > 300 GeV 


°-'-1.0 


2.9 


o 9 + 1.2 
°- z -0.9 


3.2 


2 q+0.8 
z - u -0.5 


2.1 


m > 400 GeV 


2 3 +L1 
z -°-0.6 


1.7 


2 4 +o.9 

z - 4 -0.6 


2.5 


1 c+0.6 


1.8 




e e 


e /i 


A* A* 


m > 15 GeV 


23.2lf:| 


25.7 


26.2l^° 6 6 


34.4 


12-1111 


18.5 


m > 100 GeV 


12.01U 


18.7 




16.9 


6.0l 2 j 


10.1 


m > 200 GeV 


4.91H 


4.0 


4.6l 2 i 


4.5 


2 7 +i.i 

z - ' -0.7 


4.4 


m > 300 GeV 


2 9 +L0 
z - y -0.6 


2.7 




3.5 


1 c:+0.8 
l-°-0.3 


1.7 


m > 400 GeV 




2.3 


o q+0.8 
z '°-0.5 


2.5 


-i 9+0.4 

1 - z -o.o 


1.2 



two isolated leptons with large transverse momentum within the pseudorapidity range of 
the tracking system (|r/| < 2.5). The 95% confidence level upper limits on the cross section 
of new physics processes within this fiducial region range between 1.7 fb and 64 fb for l^l^ 1 
pairs depending on the dilepton invariant mass and flavour combination. 
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